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Abstract 

The infrared behavior of the gluon propagator is directly related to confinement in 
QCD. Indeed, the Gribov-Zwanziger scenario of confinement predicts an infrared 
vanishing (transverse) gluon propagator in Landau- like gauges, implying violation 
of reflection positivity and gluon confinement. Finite-volume effects make it very 
difficult to observe (in the minimal Landau gauge) an infrared suppressed gluon 
propagator in lattice simulations of the four-dimensional case. Here we report results 
for the SU(2) gluon propagator in a gauge that interpolates between the minimal 
Landau gauge (for gauge parameter A equal to 1) and the minimal Coulomb gauge 
(corresponding to A = 0). For small values of A we find that the spatially-transverse 
gluon propagator D tr (0, \p\), considered as a function of the spatial momenta \p\, is 
clearly infrared suppressed. This result is in agreement with the Gribov-Zwanziger 
scenario and with previous numerical results in the minimal Coulomb gauge. We 
also discuss the nature of the limit A — > (complete Coulomb gauge) and its relation 
to the standard Coulomb gauge (A = 0). Our findings are corroborated by similar 
results in the three-dimensional case, where the infrared suppression is observed for 
all considered values of A. 
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1 Introduction 



The infrared behavior of the gluon propagator is linked to the confinement of 
gluons.[l] In particular, the confinement scenario of Gribov and Zwanziger[2,3,4,5,6,7,8] 
predicts a (transverse) gluon propagator vanishing at zero (Euclidean) momen- 
tum in Landau gauge and in Landau-like gauges (or A-gauges). The latter class 
refers to gauges interpolating between the Landau and the Coulomb gauge [9], 
with a gauge condition (in d dimensions) given by 

Xd A a + diA" + ... + d^A^, = , (1) 

where the gauge parameter A is between 1 and 0. Let us recall that an in- 
frared (IR) null gluon propagator has far-reaching consequences. Indeed, such 
a particle cannot have a positive semi-definite spectral function[10,ll,12,13] or, 
as a consequence, a Kallen-Lehmann representation. This is regarded as one 
possible manifestation of confinement, [1,14, 15] when considering Euclidean 
correlation functions!"*! 

The question of whether the Landau-gauge gluon propagator is indeed null at 
zero momentum is a long-standing one. Various continuum methods, based on 
functional approaches, yield a vanishing gluon propagator. [1,6,7,8,16,17,18,19,20,21,22,23,24,25] 
This result is rather tightly constrained, [25] i.e. it seems to be the only possible 
solution satisfying both Dyson-Schwinger equations and functional renormalization- 
group equations. At the same time, lattice calculations in four dimensions 
have obtained an IR-suppressed Landau-gauge gluon propagator D(p) only 
when using strongly- asymmetric lattices [26] or a coupling constant in the 
strong-coupling regime. [27] For lattice couplings in the scaling region and us- 
ing symmetric lattices, one finds for the Landau-gauge gluon propagator an 
increase slower than 1/p 2 as one approaches the IR region[28,29,30,31], with 
a finite value for p = 0. [32,33,34,35] The fact that a propagator decreasing at 
small momenta is not observed in the 4d case, even for volumes of almost (10 
fm) 4 , [35] is probably due to very strong finite-size effects. [32,33,34,35,36,37,38,39] 
This assumption is supported by numerical results in the 3d case, where much 
larger lattice sides are accessible. In this case, there is substantial evidence 
for an IR-suppressed gluon propagator in Landau gauge, [40,41,42] in agree- 
ment with continuum calculations. [6,43] However, also in this reliable 
extrapolation of -D(O) to the infinite- volume limit is still lacking. [42] Finally, 
let us recall that lattice Landau calculations[44,45,36,37] have also obtained 
direct evidence for the non-positivity of the gluon spectral function, both in 
the three- and in the four-dimensional cases. 

The Gribov-Zwanziger confinement scenario applies also to Coulomb gauge. [2,3,4,46,47,48,49,50] 
In this case it is important to observe that the standard Coulomb-gauge-fixing 

1 Thus, it is a sufficient condition for gluon confinement. [1,15] 



2 



condition d\A\ + . . . + dd-iA^_ x = is not a complete one, due to the resid- 
ual gauge degrees of freedom g(t). On the other hand, a possible complete 
Coulomb gauge condition [9] can be obtained using the class of gauges defined 
in (1). Indeed, the parameter A interpolates between the Landau (A = 1) and 
a complete Coulomb gauge, corresponding to the limitO] A — ► 0. Therefore, 
the complete Coulomb gauge condition is, by definition, a smooth limiting 
case of the interpolating gauge (1) while, of course, this is not the case for the 
standard Coulomb condition (A = 0). Let us recall that the gauge condition 
(1) above can be obtained by minimizing the (lattice) functional 3 1 



%] = -Tr£ 



d-l 

XU (x) + J2 u i(x) 

i=l 



where U^x) indicates a lattice link variable in the /i direction. Then, the lim- 
iting case A — > corresponds to minimizing the following two functionals[51] 



W^)] = - Tr E £ Ufa) (3) 

X 8=1 

£ver[s(t)] = -Tr]T U (x) . (4) 

The minimization of the first functional is equivalent to a Landau gauge 
condition fixed on each time slice, using g(x) gauge transformations, i.e. it 
corresponds to the standard (incomplete) Coulomb gauge. The minimization 
of the second functional, considering only g(t) gauge transformations, pro- 
vides additional constraints, necessary to eliminate the residual gauge degrees 
of freedom. Note that we can also write £ver[<?(£)] = — Tr J2t Qo(t) , with 
Qo(t) = J2s U (t,x) . Then, the minimization of £ver is like a one- 
dimensional Landau gauge fixing. Of course, quantities defined in terms of 
spatial link variables Ui(x) are not affected by the residual gauge condition 
obtained by minimizing the functional £ver [#(£)]• 

Numerical studies in minimal Coulomb gauge have shown [for the SU(2) case 
in 4d] that the instantaneous transverse gluon propagator D tr (p) is indeed 
suppressed in the IR limit. [51,52,53,54,55] Also, in the infinite-volume limit, it 
has been found [5 1,52, 53] that D tr {p) is well described by a Gribov-like prop- 
agator with a pair of purely imaginary poles m 2 = ±iy. These results are 
in agreement with the Gribov-Zwanziger confinement scenario. [2,3,4,46,47,48] 
The fact that, for a given lattice size L (in fm), one sees an IR-suppressed 
transverse gluon propagator in 4d-Coulomb gauge and in 3d-Landau gauge 
but not in 4d-Landau gauge may be related to a quantitatively different IR 



2 Clearly, since the gauge fixing (1) is complete for any A ^ 0, it is also a complete 
one when considering the limit A — > 0.[9] 

3 A similar functional can be defined in the continuum. 
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suppression in the two cases. Indeed, functional methods[6,43,56,57,58,59] pre- 
dict a stronger suppression in 4d-Coulomb gauge and in 3d-Landau gauge than 
in 4d-Landau gauge, i.e. the so-called IR gluon exponent should be larger 
for the 4d-Coulomb and the 3d-Landau cases. 

One should note that, for any non-zero value of A, the gauge condition (1) is 
essentially a deformed Landau gauge, i.e. the IR exponents of the propaga- 
tors do not depend on A [60]. In particular, calculations using Dyson- Schwinger 
equations suggest [60] that, at momenta sufficiently small compared to a sep- 
aration momenta p s , the transverse gluon propagator behaves as in Landau 
gauge, i.e. all Lorentz and color components of the gluon propagator vanish 
at zero four-momentum. However, the limit A — > can also be thought of as 
sending to zero the momentum p s , which separates Coulomb-gauge-like from 
Landau-gauge-like behavior. Indeed, for any finite value A ^ and given the 
Landau gauge condition d^A® = 0, we can obtain the gauge condition 

d Q A a + A" 1 [d x Al + . . . + da^Al^] = (5) 

by using the rescalinsp"! X; — > Xj/A for i ^ 0. This implies, in momentum space, 
the rescaling pt — > \pi. Thus, if we consider only spatial momenta, we find 5 1 
that p s is rescaled to Xp s and goes to zero when A — > 0. As a consequence, one 
should expect that, for very small values of A, all correlation functions would 
show a Coulomb-like behavior for momenta p > p s , with p s very small. In 
particular, the correlation function that corresponds to the transverse (instan- 
taneous) gluon propagator in Coulomb gauge should become more and more 
IR suppressed as the parameter A approaches zero. Investigating whether this 
is the case is the aim of this work. 



2 Numerical results 



Following Eqs. (8) and (9) of Ref. [60], we can consider on the lattice the 
three-dimensionally transverse gluonic correlation function 

At zero four-momentum one has 

< Af(0)A%0) > 
D tr 0) = S J 7) 



4 In Ref. [60], a similar rescaling was used to show that Dyson- Schwinger equations 
for A-gauges are equivalent to the Landau case for all A ^ . 

5 Of course, this simple explanation is correct only at tree-level and can be modified 
by the renormalization of A. 
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Here, N c is the number of colors, V is the d- dimensional lattice volume, Lorentz 
indices i,j are summed only over the d — 1 spatial directions and A%(p) is 
the d- dimensional Fourier transform of the gluon field. Let us recall that, by 
considering only spatial momenta (i.e. p = 0), the function D tr (0, \p\) is 
predicted to be IR suppressed — and vanishing at zero momentum — for all 
non-zero values of A in three and in four dimensions. [60] Also, when A is null, 
the above definition yields the instantaneous part of the three-dimensional 
transverse gluon propagator [51] 



where now the Fourier transform of the gluon field is evaluated for each time 
slice. Indeed, when A is null, the gauge transformations g(t) are independent of 
the gauge transformation g(x) and the two sets of transformations commute. 
Then, the terms in Eqs. (6) that depend explicitly on g(t) are averaged to zero 
and one is left with the expression above. Note that this explanation is valid 
whether the residual gauge freedom g(t) is fixed or not. 

Here we evaluate numerically _D tr (0, \p\) as a function of \p\ for SU(2) Yang- 
Mills theory in three and in four dimensions for several values of the parameter 
A. Details of the simulations can be found in. [42,61] Let us note that for 
A ^ 1 the numerical gauge fixing is very similar to the usual Landau gauge 
fixing. [41] On the other hand, when A goes to zero one sees[62] that more 
iterations are needed in order to satisfy a given numerical accuracy for the 
gauge fixing, in agreement with a recent analytic study. [63] This problem can 
be partially reduced by adding some extra gauge fixing sweeps in the \i = 
direction for each iteration of the gauge-fixing algorithm, i.e. by considering 
gauge transformations g(t) that depend only on the fi = component of x. 
Finally, we did not consider here possible systematic effects related to the 
breaking of rotational symmetry or to the existence of Gribov copies. The 
former type of effects can be parameterized by [64] a 2 p^\ where a is the lattice 
spacing and = yj^pfj. Therefore, this type of effects are not expected 
to play a significant role in the IR region, considered here. As for the latter 
type of effects, in the infinite-volume limit, averages taken over configurations 
belonging to the so-called Gribov region Q should coincide [8] with averages 
obtained by restricting the functional integral to the so-called fundamental 
modular region T, whose interior is free of Gribov copies. 

Our results^ are reported in Figure 1. In three dimensions, a well-defined 
maximum (and thus an IR suppression) is visible for all values of A. This 
includes also Landau gauge (A = 1), confirming earlier results. [40,41,42] As 
can be seen from the plot, the maximum value attained by the propagator 
seems to move to larger momenta with decreasing A, going from about 400 

6 Preliminary results have been presented in Ref. [65,66]. 
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D tr as a function of |p| in three dimensions 
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D tr as a function of |p| in four dimensions 
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Fig. 1. The gluonic correlation function D tr (0, \p\) in three (top figure) and four 
(bottom figure) dimensions. The spatial momenta \p\ are chosen along a single axis. 
Circles indicate data for A = 1 (Landau gauge), crosses are used for A = 1/2, squares 
for A = 1/10, triangles for A = 1/20, stars for A = 1/100 and upside-down triangles 
represent results at A = (Coulomb gauge). Data have been obtained at (3 = 4.2 in 
three dimensions and at (3 = 2.2 in four dimensions. The lattice size is 40 3 ~ (6.9 
fm) 3 in three dimensions and 22 4 ~ (4.6 fm) 4 in four dimensions. In addition, in 
the bottom panel, diamonds correspond to a 40 4 ss (8.4 fm) 4 lattice at A = 1/100. 
The physical scale has been set using Refs. [42,67]. Note that our error bars are 
smaller than the sizes of the symbols. Also, for A = 0, we consider the standard 
(incomplete) Coulomb gauge. 
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D 00 as a function of |p| in three dimensions 
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D 00 as a function of |p| in four dimensions 




Fig. 2. Same as in Fig. 1, except that here we show data for the Ajo(0, \p\) tensor 
component. 

MeV in the Landau case [42] to about 600 MeV for the Coulomb case. At 
the same time, as A decreases, the maximum becomes more visible. In four 
dimensions, on the other hand, no discernible peak is visible in Landau gauge 
when considering the volume V = 22 4 . Decreasing A, however, leads to a 
suppression of the propagator in the IR region. In particular, at the smallest 
value of A considered, i.e. A = 1/100 and for a lattice volume V = 40 4 , a 
maximum is seen also in four dimensions. This maximum is as visible as the 
one in three dimensions in Landau gauge. Thus for small A, just as for Coulomb 
gauge, [51] one sees a maximum of the transverse gluon propagator already for 
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relatively small volumes. 



A similar A- dependence is also seen in the tensor component Dqq(0, \p\). In 
particular, as A becomes small, this tensor component of the gluon propaga- 
tor becomes enhanced for small (nonzero) values of \p\ (see Fig. 2). Let us 
recall that, in Coulomb gauge, the D Q0 component is IR divergent already 
at the perturbative level. Moreover, the instantaneous part of Dqq is related 
to the color Coulomb potential and should be diverging as l/|p| 4 at small 
momenta. [46,47,48,51,52,53,54,55] Thus, from the discussion presented in the 
Introduction, one could expect to observe a D 00 component enhanced at inter- 
mediate momenta p > p s but still suppressed at small momenta. This is indeed 
the case and from the plots in Fig. 2 one clearly sees how this enhancement 
is developing for \p\ ^ as the value of A decreases. Let us stress that the 
discontinuity in the behavior of the propagator at \p\ ^ is observed since we 
are taking the limit A — > 0. A similar discontinuity is also obtained when im- 
plementing the complete Coulomb gauge described in the Introduction^ On 
the other hand, for the incomplete Coulomb gauge, D 00 is clearly a continuous 
function of \p\ (see Fig. 2). 

As said in the Introduction, the IR exponent for the transverse correlator 
is predicted to show a discontinuity in the limit A — > 0. This is not seen from 
our data. However, a reliable check of this prediction can only be obtained 
if one has control over the infinite-volume and the continuum limits, which 
we have not yet achieved. On the other hand, taking the limit A — - > clearly 
induces a discontinuity in the behavior of the temporal component .Doc [60] 



3 Conclusions 

We have presented the first direct observation of an IR-suppressed gluonic 
correlation function on a symmetric 4d lattice in a Landau-like gauge, with 
gauge parameter A as in Eq. (1). (Landau gauge is obtained for A = 1, while 
Coulomb gauge corresponds to A = 0.) The suppression is seen for sufficiently 
small A when considering moderately small lattice volumes. Judging from the 
results shown for the 3d case, it is conceivable that a similar suppression might 
be observed for any A if a large enough lattice side is considered. Furthermore, 
since the limit A — > 1 is smooth, [9,60] we expect to see an IR suppression for 
sufficiently large lattices also in Landau gauge. Moreover, we have obtained 
that the Dqq component of the gluon propagator gets enhanced in the IR region 
(for \p\ ^ 0) at small values of the interpolating parameter A. These results 
clearly demonstrate the transformation of the correlation functions from Lan- 
dau gauge towards Coulomb gauge when considering the limit A — > 0. They 

7 This is a consequence of the minimization of £ V er [#(£)] , defined in Eq. (4). 
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also provide additional support to the Gribov-Zwanziger scenario of confine- 
ment, establishing an IR suppression of gluonic correlation functions in four 
dimensions beyond Coulomb gauge. The use of the interpolating gauge thus 
constitutes a promising alternative to studies in Landau gauge, since finite- 
size effects are significantly smaller, allowing a more efficient investigation of 
the gauge dependence of correlation functions of confined objects. 
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